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In the present work we study the two-point correlation function R{e) of the quantum mechanical 
spectrum of a classically chaotic system. Recently this quantity has been computed for chaotic and 
for disordered systems using periodic orbit theory and field theory. In this work we present an in- 
dependent derivation, which is based on periodic orbit theory. The main ingredient in our approach 
is the use of the spectral zeta function and its autocorrelation function C(e). The relation between 
R{e) and C(e) is constructed by making use of probabilistic reasoning similar to that which has 
been used for the derivation of the Hardy - Littlewood conjecture. We then convert the symmetry 
properties of the function C(e) into relations between the so-called diagonal and the off-diagonal 
parts of R{e). Our results are valid for generic systems with broken time reversal symmetry, and 
with non- commensurable periods of the periodic orbits. 



PACS numbers; 03.65.Sq, 05.45. -|-b 



I. INTRODUCTION 

Quantum chaologyS has attracted thft,attention of the 
physics community after the discoverycl that the spec- 
tral correlations of classically chaotic systems are univer- 
sal. They exhibit a strong level repulsion, which induces 
a non-trivial two-point correlation function R{e) of the 
density of states. 

The spectral rigidity of a generic chaotic system was 
computed for the first time by BerryB. He expanded the 
density of states over periodic orbits by making use of 
Gutzwiller's trace formulau. He estimated number of 
orbits of a given length from the Hannay - Ozorio De 
Almeida sum ruleB and computed the so-called diagonal 
part of the form-factor i^TdiagCr). The form-factor K{t) 
is a function of time, and is the Fourier transform of the 
two-point correlation function R(e). 

The semiclassical theory of the form-factor distin- 
guishes between three main time scales, see diagram 
Fig. |l|(a). The contribution of the short periodic orbits 
is shown schematically as a sequence of (5-peaks, and the 
relevant time scale is the period of the shortest periodic 
orbit Tp. This non-universal behavior of the forra.- factor 
prevails between r = and the ergodic time Toig.tl Berry 
found that K{t) oc t for tr t ^ ''org, and showed that 
this explains the level repulsion. 

In the present work we are primarily interested in the 
semiclassical (periodic-orbit) theory of the form-factor 
for T ^ Th, where th is the Heisenberg time. This is 
the third time scale. It is of quantum mechanical nature, 
because the Heisenberg time is proportional to the mean 
density of states. In Fig. |^ we have chosen units of time 
such that Th = I. Due to the discrete nature of the den- 
sity of states, the form-factor becomes constant for times 



much larger than th, see Fig. 0(a). The semiclassical 
theory for the form-factor for t ^ th iias been analyzed 
recently by Bogomolny and Keating.Ll They found the 
fingerprints of the short periodic orbits in the vicinity 
of the Heisenberg time as shown in Fig. [|(a) by (5-peaks 
near t = th- 

The most difficult problem addressed by Bogomolny 
and Keating was the computation of the off-diagonal part 
of tb£ two-point correlation function i?off(e)- Argaman 
et aE pointed out that it is possible to compute i?off(e) if 
the two-point statistics of the classical actions is known. 
Let us call the set of actions of all primitive periodic 
orbits (PPOs) the length spectrum. The density of states 
is related to the length spectrum by the Gutzwiller trace 
formulaU. Therefore correlations of eigenenergies should 
be related to correlations of actions. 

Our approach to the semiclassical evaluation of fioff (s) 
further develQps-.the approach of Argaman et aU and 
their followeraJ'liJ but makes use of a different starting 
point. The first new element is the study of the spectral 
zeta function (^{E), its autocorrelation function C{e), and 
its Fourier transform K'^{t). 

The spectral zeta function C,{E) is an important tool 
for an analytic or numeric computation of the energy 
levels of any quantum mechanical system. The zeros of 
this function are the eigenenergies of the HamiltoniaxL. 
A unique definition of C,{E) is to be given later. VoroslliJ 
has proposed to compute this function by making use pL 
a product over the periodic orbits. Berry and KeatingEj 
have expanded this product over the composite periodic 
orbits (CPOs), see precise formula below. The spectral 
zeta function is a smooth function that has no (5-peaks 
like the density of states does. Therefore it is a good idea 
to characterize level statistics of chaotic systems by the 
C,{E) autocorrelation function. This was done by Kette- 
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FIG. 1. Form-factors derived from correlation functions of density of states and spectral determinant. Diagrams (a), and (c) 
correspond to a small chaotic system, whore one can observe fingerprints of short periodic orbits in the form-factors, shown as 
(5-functional peaks. Diagrams (b) and (d) correspond to diffusive systems, which can be considered as huge complex chaotic 
systems. One can observe smoothing of the form- factors near r = and t = th- Particular shape of this smoothing contains 
system specific information like dimensionality and characteristics of disorder. 
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mann, Klakow and Smilansky, who also computed this 
autocorrelation function C{e) for two and three dimen- 
sional Sinai billiards .113 

The diagram Fig. 0(c) schematically shows K''{t). For 
systems with broken time reversal symmetry, as in Fig. ^, 
one can see the separation of the time scales very clearly. 
The short time behavior of K''{t) is determined by the 
short CPOs, and it is represented by (5-peaks. They 
become dense giving a constant behavior of K''{t) for 
T > 7'org, see Ref. |l3[ It is known that the spectral 
zcta function has to satisfy a functional equation.EJ This 
equation can be derived from the definition of C(^) ^nd 
implies an exact mirror symmetry of K''{t) around the 
half Heisenberg time: K'^{th/2 + t) = K'^{th/2 - r). 
For this reason one has to observe the fingerprints of the 
short composite periodic orbits near th, they are shown 
as 5-peaks near th in Fig. |l|(c). Therefore, if the short 
time behavior of K'^{t) is known, then the behavior of 
K'^ir) near th is also known. 

Similar to the case of R{e), the periodic orbits ex- 
pansion of C(^) allows one to separate C{£) and 
into diagonal and off-diagonal parts. Following Berr; 
we can assume that the off-diagonal part of K^(t) van- 
ishes for T < Tjy/2. This assumption together with the 
mirror symmetry of K'^{t) implies an explicit connec- 
tion between the diagonal and the off-diagonal parts of 
K''{t). Physically, this relation has to be interpreted as 
a quantum- classical time scale separation. 

Since C(-^) is expressed in terms of CPOs by the Berry 
- Keating formulallj, the autocorrelation function C(e) is 
related to correlations in the composite length spectrum., 
which is the set of all CPOs. In this way we introduce an 
analog of length correlations, which were used by Arga- 
man et aB to discuss the two-point spectral functions. 
One of the central ideas of the present work is to com- 
pute the correlations in the composite length spectrum by 
making use of the mirror symmetry of K''{t). More over, 
the time scale separation provides an explicit dependence 
of K''{t) on Th- The two point correlation hmction of 
composite actions is just the Fourier transform of K'^{t) 
with respect to hidden in th- This idea is similac to 
one which has been suggested by Balian and BlochJlj 

The length or the action of the given composite pe- 
riodic orbit is simply the algebraic sum of the lengths 
or the actions of the PPOs forming this composite or- 
bit. For this reason, action correlations of CPOs have to 
be related to action correlations of PPOs. In the present 
work we construct an integral relation between these cor- 
relation functions by making use of simple probabilistic 
arguments. This relation is generic and independent of 
the dimensionality of the system, however it might be 
dependent on the symmetry of the system. We perform 
calculations only for systems which does not have spa- 
tial symmetries, because we assume that the composite 
length spectrum is non-degenerate. 

The statistical relation between correlations of CPOs 
and correlations of PPOs can be converted by the inverse 
Fourier transform and a special regularization procedure 



to a relation between the off-diagonal parts of K{t) and 
K''{t)- Calculations show that the behavior oi K{t) near 
Th, which is shown schematically in Fig. Ma), reproduces 
the behavior of K'^ir) near th, see Fig.^(c). 

The ideas of this paper are best understood by the 
following logical fiow: 

1. The off-diagonal part of K{t) is the Fourier trans- 
form with respect to \/h of the two-point statistics 
of actions of PPOs. 

2. The two-point statistics of actions of PPOs can be 
computed from the two-point statistics of actions 
of CPOs. 

3. The two-point statistics of actions of CPOs is de- 
termined by the off-diagonal part of i^^(r). 

4. The diagonal and off-diagonal parts of K'^{t) are 
connected, because C(^) satisfies the functional 
equation. 

5. For the systems with broken time reversal symme- 
try there is a clear separation between the classical 
scale (time of mixing) and the quantum mechanical 
scale (time of quantum recurrence) . In this case we 
can compute the off-diagonal part of K"^ (r) explic- 
itly, if the diagonal part is known. 

6. The diagonal parts of either K{t) and K^ir) or 
i?(e) and C(£) have well known periodic orbit ex- 
pansions, which can be evaluated, both numerically 
and by making use of sum rules B'lJ 

The two-point statistics of actions is a property of 
classical dynamics and it is independent of fi- Unfortu- 
nately, there are only a few examples of chaotic systems 
where a sufficient number of periodic orbits is known and 
the action correlations caxi be computed directly from 
the classical Hamiltonian.EHl We circumvent this difficulty 
by computing action correlations from another source, 
namely from the diagoaal part of C(£), which is also a 
purely classical object.t^l 

The integral relation connecting the correlations of 
the PPOs and the CPOs can be applied directly to de- 
rive the well known correlations between prime numbers. 
As we show in Appendix B, the Hardy - Littlewood 
expressionllEI can be reproduced in this framework. We 
emphasize this application as a very stringent test for our 
probabilistic method. 

Large chaotic systems, for example billiards with a 
large number of scatterers or disordered systems, usu- 
ally have Terg ^ Tp, since Tcrg is the time of a diffusion 
through the system. The system specific features of K{t) 
and K'^{t) are smooth, see Figs. |l(b) and 0(d). The con- 
dition Torg ^ Tp allows us to express the diagonal part of 
C(e) as the Fredholm determinant of the diffusion propa- 
gator. In this way our expressions for R{s) reproduce prfi^ 
cisely the field theory result of Andreev and Altshuler.EZi 
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At this point it is appropriate to review the impor- 
tant work on the level statistics of disordered systems 
which had a very important impact on the develop- 
ment of t|hjG-|Scmiclassical theories presented here and 
elsewhere. 0113 Field theory was used by Efetov to com- 
pute the level-atatistics of small metallic samples, see the 
review papeill3. He considered the electron moving in 
a random potential of impurities and confined by sam- 
ple boundaries. He assumed that the fluctuations of the 
fields are uniform across the sample and obtained univer- 
sal results for the level statistics. 

Altshuler and Shklovskiill3 made use of perturbation 
theory and expressed the so-called perturbative part of 
R{e) in terms of the density-density correlation function, 
which is the propagator of the diffusion equation. Arga- 
man et ac2l showed that the perturbative part of R{e) is 
nothing but the diagonal part of R{e) mentioned before. 

The diagonal (or perturbative) part of i?(e) is singular 
at small energies; the only way to remove this singularity 
is to compute the off-diagonal (pt non-perturbative) part 
of R{s). Andreev and AltshulertZl computed this term by 
making use of non-perturbative field theory and obtained 
the answer in terms of the same density-density correla- 
tion function. These authors repeated Efetov's calcula- 
tions, but allowed the spatial fluctuations of the super- 
fields. 

The free energy functional in Efetov's theory describes 
the diffusion modes. These modes have to satisfy the dif- 
fusion equation. This theory was remarkably generalized 
by Muzykantskii and KhmelnitskiiEil. They obtained the 
free energy functional describing the eigenmodes of the 
kinetic equation. Such solutions to the kinetic equation 
have a meaning of the density-density correlationsES. 

Muzykantskii and Khmelnitskii also suggested that the 
modes of the kinetic equation should be replaced by the 
modes of the Liouvillian operator tfiOne goes from difftu. 
sive to chaotic systems. Agamei a£3 and Andreev et a£fl 
constructed field theory for chaotic systems by averaging 
over the energy and obtained the Liouvillian operator 
in the kinetic part of the free energy functional. Both 
types of field theory show that the spectral statistics of 
a chaotic system is described by the determinant of the 
Liouvillian operator. 

The connection between the diagonal part of the cor- 
relation function and the determinant of the Liouvillian 
operator is almost trivial in the framework of periodic 
orbit theory. This connection is different from the field 
theory prediction due to the terms containing repetitions 
of the PPOs. The periodic orbit theory expression for the 
off-diagonal part of the correlation function obtained by 
Bogomolny aad KeatingEl is also different from the field 
theory resulto. Therefore, the derivation of the Andreev 
- Altshuler result in terms of the action correlations is in- 
teresting, and it also gives additional information about 
the autocorrelation function of the spectral determinant. 
The lattp^ has been obtained by using random matrix 
theoryE3n3 and random polynomial theorjtZl, as opposed 
to field theory. 



All correlation lunctions employed m t ne p resent work 
are defined in Sees. || and HI. In Sec. |IV| very simple 
probabilistic arguments will help us to build the integral 
equation connecting the correlations of CPOs with the 
correlations of PPOs. The behavior of form-factors near 
the Heisenberg time is considered in Sec. 0. In the same 
section our results are compared with the universal ran- 
dom matrix theory predictions and the results of Refs. 
p^ , ^ We discuss the physics of the action correlations 
in Sec. 0. We also summarize our results in Sec. No. 



II. DEFINITIONS OF OBJECTS RELATED TO 
THE SPECTRAL CORRELATIONS. 

We start to build our theory for the specific example of 
the chaotic billiard. Let us also put the Aharonov - Bohm 
fiux through the billiard in order to break the time re- 
versal symmetry. Let us also assume that the system has 
no spatial symmetries. The classical motion of a charged 
particle in the billiard is finite and therefore one can as- 
sociate two sets with this system. TheJirst set, {Ip}, is 
formed by the lengths of all the PPOsEj, and each orbit 
is labeled by the index p. Each orbit makes fip windings 
around the Aharonov - Bohm flux. We distinguish the 
orbits with positive and negative winding numbers, i.e. 
the paths going clockwise and counterclockwise around 
the flux tube. Therefore, the length spectrum {Ip} is de- 
generate, the lengths of the orbits with non-zero winding 
numbers appear twice. 

The second set associated with the billiard is the quan- 
tum mechanical spectrum. It is formed by the wave 
vector magnitudes {fc„}. Each of them corresponds to 
the eigenvalue of the Hamiltonian £"„ = E{kn) and 
E{k) = n^k'^/{2m), where TO is the mass of a particle. 

The density of states is a sequence of 5-peaks, and it is 
usually computed by making use of the Green function 
of the system: 



d{k) = ^5{k- k„ 



1 



ImTrG"'(fc)£;'(fc) . (la) 



n=l 



Here &'{k) is the retarded Green function of the system, 
and the prime means derivative with respect to k. The 
trace of G^{k) can be represented as a sum over the perk 
odic orbits. Such an expansion is called a trace formula.cl 
The derivation of this expansion for a billiard can be 
found in the review Ref. E8. The result is 



TrGW(fc)=-EESl 



p r—1 



ind^k) , 



(lb) 



where the index p runs over the PPOs, the length of each 
orbit is Ip, the Lyapunov exponent of each orbit is Xplp, 
and the Maslov index is z/p. We also broke the time re- 
versal symmetry by adding the phase due to Aharonov - 
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FIG. 2. Schematic diagram demonstrating how to compute 
the correlation function of energy levels. Each point repre- 
sents pair of the energy levels (fc„, k„'). The correlation func- 
tion R{e, k) picks up the points inside one of the ellipses shown 
in the diagram. Three ellipses near the diagonal are drawn 
for k equal to fc'"', fe*''', and e < Afc. We draw the 

fourth ellipse for the case e Afc and the choice of relevant 
fc becomes ambiguous. 



Bohm flux (j) multiplied by the winding number of the tra- 
jectory /ip. Lyapunov exponents in Eq. (lb) are defined 
per unit length of a trajectory and not per number of 
scatterings as in Ref. This formula also contains the 
smooth part of the density of states d{k), which is pro- 
portional to the volume of the system. Equations jla] ) 
and (lb) show that the density of states contains an os- 
cillating part. It is defined as 



rfosc(fc) = d{k) — d{k) 



(Ic) 



and is proportional to the imaginary part of the first term 
on the right hand side of Eq. (lb). 

The correlation function of the quantum mechanical 
spectrum and its Fourier transform have to be defined 
in a special way, because of the (5-functional form of the 
density of states and the discrete form of the expansion 
Eq. (|rb|) . We define the correlation function by making 
use of the averaging over k with Gaussian weight 



R{e,k) 



dq 



X dosc{q + -)dosc{q - 77) 



27rAfc 
2' 



(2) 



and this definition of the correlation function is valid for 
e <C fc. We emphasize that e in our definition is the 
difference of the wave numbers and not of the energies as 
accepted. 

The definition of the correlation function Eq. (||) con- 
tains the product of two (5-functions from two densities of 
states and only one averaging. Practically, the (5-function 



remaining after the integration in Eq. (g) has to be re- 
placed by a smoothed function of the width smaller than 
the mean level spacing. Then, the definition Eq. (H) 
can be understood geometrically. The two-point corre- 
lation function is proportional to the difference between 
the number of the level pairs (fc„, fc„') inside the region, 
which we showed schematically as an ellipse in Fig. ||, 
and the square of the mean density of states multiplied 
by the area of this region. The "length" of the ellipse 
is of the order of Afc and the "width" has to be smaller 
than the mean level spacing. 

We show three ellipses demonstrating the correlation 
function computed for k equal to fc^°\ k'^^\ k^'^\ and e <C 
Afc. It is important to note that Afc has to be smaller 
than the characteristic scale of variation of i?(e, fc) as a 
function of fc. We will see later that this scale is ^ 
sd'{k), where the prime means derivative with respect to 
fc. Therefore, the choice of the averaging interval length 
is limited by the inequality ed'{k)Ak < 1. 

The fourth ellipse in the diagram Fig. |^ is drawn in 
the area where e ~ Afc. It is difficult to decide whether 
the value of fc in this case is fc^°\ or fc^**), or maybe fc*^'^^. 
In order to make the correlation function well defined for 
all e we multiplied it by Gaussian prefactor e~^^ ^ see 
Eq. d) where AZAfc > 1. 

The domain of e, where the correlation function is 
meaningful, is limited now by two competing conditions: 
elSl < 1 and e < l/{Akd'). The latter inequality has to 
be implied by the former one, and therefore we obtain 
the important condition 



Afc 



1 



< 

A/ ^ d'{k) 



(3) 



The averaging has to be performed over a large number 
of the energy levels. This number is 



AA/" - Afcd(fc) > 1 



(4) 



The inequalities Eqs. (^ and (^ are not very restric- 
tive and most of the published experimental and numer- 
ical work, where the correlation function was computed, 
employed the averaging intervals of the width satisfying 
them. 

Similar statistics can be defined for the length of the 
PPOs entering the right hand side of Eq. (|l^). The 
length spectrum of the system has the density 5{x — 

Ip), and the weighted mean density dppo(a;) is the sum 
over n,u, jj, of 

J ( \ ' A A A ^ 2A(a 

^PPO V*^5 ^5 / ^ ^ ^71 Up Up^ fl flp ' 



P 



2ttAI 



rfppo(-'^) , X ^ {x, n, V, ^) 



(5) 



where we have set the averaging interval to be precisely 
A^, because we are going to build analytical relations 
between the statistics of energy levels and the statistics 



5 (November 30, 1997) 




Ja) — 



Ip 

FIG. 3. Schematic diagram demonstrating how to compute 
the correlation function of the length spectrum. Each point 
represents the pair of the primitive periodic orbits lengths 
{lp,lpi). The correlation function R{x,y) picks up the points 
inside one of the ellipses shown in the diagram. The three 
ellipses near the diagonal were drawn for x equal to x^°'\ x^''\ 
x^'^\ and y <^ Al. The fourth ellipse was drawn for the case 
y Al and the choice of relevant x becomes ambiguous. 



of periodic orbits. The mean density in Eq. picks 
up the orbits p with the defined value of the number of 
the wall reflections Up. We keep this number in Eq. (|^), 
because the action correlations were found between the 
orbits with the same Up, see Ref. 

The correlation function of the length spectrum is de- 
flned in a slightly different way from the correlation func- 
tion of the density of states: 

R2ix,n,iy,^j.;x',n,iy\fj,') = Snn^S^ ^,^6^ 

p 



p' 



2ttAI 



^ ^ e ^ ^n' n^/^u' u^i^^' ^^/^{x X ~^ ^p') 
p'^p 

I x + x' 'p + 'p' ^2 

rv 2 11 



e 2AH 



2-kAI 
R2{X;X') , 



(6a) 



where the first term in the braces is the product of the 
mean densities of PPOs and the second term is the prob- 
ability of finding two PPOs having the length difference 
y. Both terms are weighted by stability factors which 
compensate the exponential proliferation of the PPOs. 

The key assumption of the present theory is that the 
orbits with different Maslov indexes or different winding 
numbers or different number of scatterings do not con- 



tribute to the correlation function of the length spectrum, 
which we should define as 

R{x, y) ^ ^ R2ix+^,n,i^,^i;x-^,n,i^,n) . (6b) 

Therefore we assume that R2{X; X') decay very fast with 
|n — n'l, — j/'l, and |/i — fi'\. The correlation function 
R2{X; X') is meaningful for |a;— ^ Al and this is taken 
into account by the factor g-^'^^i^-^'f/'^ in Eq. (|6^), 
since 1/Ak < as we will show later. 

The correlation function defined by Eq. ( |6b| ) has a nor- 
malization condition which is an extremely useful tool for 
verifying results. We have from Eqs. (^) 

J2 J dx' R2iX;X')=dppoiX) (7a) 
dy R{x,y) = dppoix) (7b) 



which are valid under the condition Al <^ x. In order 
to derive this relation we have replaced J2p^p'^p 
Eq. ( |6a| ) by two terms J2pp' ~ ^pp' ■ The second term 
gives the mean density on the right hand side of Eqs. (|^). 
The first term is a product of the mean densities and can- 
cels another term of the type ^q. (|6^. 

The correlation functions defined by Eqs. m) have the 
graphical representation in Fig. ^ quite similar to that in 
Fig. ||. In the case of Fig. ^ we are interested in the num- 
ber of PPO pairs whose mean length is near x and dif- 
ference of lengths is near y. Correlation function counts 
the number of PPO pairs inside one of the ellipses shown 
in Fig. ||. The "length" of the ellipses is determined by 
the width of the Gaussians in Eq. ( |6a| ) and to be ~ AL 
The "width" of the ellipses does not appear explicitly in 
Eq. (pa]), meaning that S{y — lp + lpi) has to be smoothed 
to the order of l/dppo(a;). 

The length of the averaging interval, Al, has to be as 
large as possible, but smaller than the characteristic scale 
of the variation of R{x, y) as a function of x for constant 
y. This scale can be estimated from the results of Ref. ^ 
and it is d'{k)/y taken for k = kx, where kx is the inverse 
density of states: 



X = 2nd{kx) , 



(8) 



which is a classical quantity. The choice of the av- 
eraging interval is, therefore, limited by the inequality 
A^ < d'{k)/y. This inequality has to be fulfilled in the 
interval of y, where the correlation function Eq. (^) is 
meaningful, i.e. for < y < 1/Afc. Therefore, the choice 
of the averaging window width is limited by the condition 
analogous to Eq. (0) 



(9) 



This inequality holds regardless of Eq. (|^) because one 
can study the statistics of energy levels independently 
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from the statistics of periodic orbits. From this point of 
view the inequality Eq. (|p) is not very restrictive either. 

The statistical properties of the PPOs length spectrum 
arc ultimately connected with the st|atistical properties 
of the quantum mechanical spectrumlj. We obtain from 
Eqs. (|l]) the periodic orbit expansion of the spectral cor- 
relations form-factor 

r°° de 

J-oo 



47r2 



pr 



e 2A1- 



(10a) 



e 2 



Kos{x, k) 



p V _ 1 r I T I ■ 
i e 2aI^>-i 



(10b) 



47r2 



^ g»fc(;p-;p,)-Afc2(ip-ip,)V2 



X e 



*(i'p-i^p')'r/2+i0(Mp-A'p') 



-^x" I dye^'^yR{\x\,y) 



(10c) 
(lOd) 



which is valid for k ^ Afc. Under this condition wc 
have neglected the Fourier transform of the first term in 
the braces on the left hand side of Eq. (|6a|), because it 
contains the factor e"'' /(^Afe )^ 



Equation (10c) is justified for x such that 



X\x\ > 1 



(11) 



where A is the mean Lyapunov exponent defined per unit 
length of the trajectory. This condition allows us to ex- 
pand the denominators in Eq. (^) and to neglect the 
summation over the repetition index r. We also assumed 
that the magnetic flux (j> is so large that (fiAfj, ^ 1, where 
Afi is the characteristic scale of /i — /^'-dependence of 
R2iX;X'). 

The key question of the theory is how to perform the 
summation over the periodic orbits if thalength spectrum 
is unknown. rThis question was resolvedu with the help of 
the sum ruleO, valid for ergodic systems. One can deduce 
the density of PPOs from this sum rule: 



^ppo (^) 



(12) 



This expression is valid asymptotically for x larger than 
the ergodic length (rorg times velocity) and it gives 
Kdiiigix) cx x. 

The definition of the form-factor, Eqs. (p^, is not com- 
plete, because we did not specify yet how to choose the 
sizes of the averaging windows Afc and Al. Common 
wisdom was that Eqs. (^0|) would survive if they do not 
break the conditions 



Ale^''/x > 1 , Afcd(fc) > 1 



(13) 



We want to define the form-factor in such a way that 
it will carry information about the correlations of en- 
ergy levels and about the correlations of periodic orbits. 
Therefore, we require Afc and Al to satisfy Eqs. (^ and 
(0) i.e. 



d'{k) < ^ < d'(fc.) . 



(14) 



Near the Heisenberg length x ^ 27rc?(fc), we have k^ = k 
and two inequalities in Eq. (|l4|) give us AZ/Afc ~ d'{k). 



III. DEFINITIONS OF OBJECTS RELATED TO 
THE SPECTRAL ZETA FUNCTION. 

One can construct an infinitely large number of com- 
plex analytic functions of fc having zeros only at fc = fc„ . 
Each such function can be considered as a spectral deter- 
minant. We are going to define the spectral determinant 
in terms of the retarded Green function G"'(fc). The trace 
of this operator is an analytic function and we have 



C(fc) = e/o'°''''{'^''^''('')^'(''^+"'^('^)} 



(15) 



This definition implies the functional equation for the 
spectral zeta function 



C(fc) = e'"r^«'^"(«'C*(A:) 



(16) 



where the asterisk stands for the complex conjugation. 
The expression in braces in Eq. ( [Tsl ) can be written as 
the series expansion over the periodic orbits Eq. ([T^). 
The epqponential function of this series was computed by 
VorosM and the result is 



C(fc) 



nn 

p r— 



^i{klp-i^p-iT/2+(pfip)^-{r-\-l/2)Xplp 



(17) 



This expression is taken from Eq. (6.21) of Ref. and 
Xp here is the Lyapunov exponent per unit length. Here 
we added the phase due to the Aharonov - Bohm flux. 

The autocorrelation function of the spectral zeta is de- 
fined similarly to the correlation function of the density 
of states: 



C(e,fc) 



V2 



dq 



/27rAfc 



xC('z + f)C(?-|) 



(18) 



Equation. (^6|) is also valid for the correlation function 

C{e, fc) = e2'^*'^"('=)'^C(-e, fc) , (19) 

where we expanded the integral in the exponent in 
Eq. ( p^ under the condition e <C fc, and also used the 
inequality Eq. (0). 
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As usual, we put the prefactor in Eq. (|18|), which will 
be translated into the averaging over the lengths of the 
CPOs. The CPOs appear from the evaluation of the 
product in Eq. and we define them as all possible 
sets of the PPOs taken without repetitions: 



c = {p} and Ic — Ip 



(20) 



The CPOs have the composite Lyapunov exponents 
Ac = T^pec'^p^p^ composite Maslov indexes Vc ^ 
Epsc'^P' composite number of wall reflections Uc = 
Spgc^P' the number of the primitive components uric = 
Spec-*^' and the composite winding numbers /j,c = 
^pGcf^p- This definition of the CPOs leads to the high 
degeneracy of the composite length spectrum {Ic}- 

The composite length spectrum has the density 
J2c^i^ ~ 'c) and the weighted mean density dcpo{x) is 
the sum over n,v, n of 



c'cpo(^) 



(21) 



The correlation function of the composite length spec- 
trum is defined similarly to Eq. (^) 

C 



Tic ^i^' i^f^t ^fi' fi^/ ' 



2ttAI 



/27rA? 



> , e 5n' n^, 5u' v^, 6{X - x' - lc + Ic' ) 



-(- 



2ttAI 



C2{X-X') , 



(22) 



^"(2;, y) = E <^2(a; + |, n, I/, ^; a; - -, n, i^, /i) . (23) 

We will assume in what follows that C2{X\ X') falls off 
rapidly with |n — and — The correlation 
function C(x,y) is meaningful only for y ^ AZ, and this 
is taken into account by the Gaussian prefactor. The 
normalization ofthe correlation function can be deduced 



from Eqs. 1^ and (g2| 



E 

n' ly' fi' 



dx'C2{X;X')^dcpoiX) , 



dyC{x,y) = dcpo(x) , 



(24a) 



(24b) 



and it is valid under the condition Al <ti x. The deriva- 
tion is similar to the derivation of Eq. (0). 



The statistical properties of the CPOs length spectrum 
are connected with the statistical properties of the spec- 
tral zeta function. We can see that by expanding the 
Fourier transform of the correlation function Eq. (|T^ ) 
over the CPOs 

-e--C(£, k) = Ki^^ix) + kU^, k) 
-00 ^~ 

(25a) 

^diag(2^) = E E i n n'^^id -e-^p'pJ)2 ( 



2ttAP 

^cpo (^) 



(25b) 
(25c) 



27rA/2 



X gifc(;c-ie')-Afe^(ic-ie')V2gi7r(mc-me) 



dye'^^yCix^y) 



(25d) 
(25e) 



where the repetiti on in d exes rp in Eg. ( ^5b| ) run from 1 
to 00. Equations (ESq), (25d) and (25c) are valid under 
the condition Eq. (|1 1|) . This condition allows us to keep 
only the factors with r = in Eq. ([l7|), and therefore 
our definition of the CPOs Eq. ( |20| ) does not contain 
the repetitions of the PPOs. We also assumed that the 
magnetic flux (f> is so large that 0A/i 3> 1, where A/i 
is the characteristic s cale of the /i — /i'-dependence of 
C2{X; X'). Equation (25a) deflnes the form-factor ofthe 
autocorrelation function of the spectral zeta. It is zero 
for negative x, whereas K(x) depends on |a;|. 

The avemged density ofthe CPOs length spectrum was 
computedtJ for a family of strongly chaotic systems and 
its asymptotic form in the general case is 



4po(a:) ^ 7 



(26) 



where 7 is the normalization constant, and 7"^ is of the 
order of the billiard size. This density is sufficient for 
the c omp utation of the diagonal part of the form-factor 
Eq. (25c) but the non-diagonal part can be computed 
only if the correlations in the CPOs are known. 



IV. PROBABILISTIC THEORY OF THE 
LENGTH SPECTRUM CORRELATIONS. 

The length spectrum of PPOs is a part of the length 
spectrum of the CPOs. Therefore we can ask what is 
the probability to find a PPO in some interval of the 
composite length spectrum. Let P be the set of all PPOs 
and let C be the set of all CPOs. Then P C C. It 
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may happen that c G C consists of only one PPO and 
therefore c € P. We define the probabiHty of this event 



Pic) 



ceP 



(27) 



It is a function of Xc where Xc = {lc,nc,Vc, y^c)- The 
probabihty P{c) counts CPOs, which are actually PPOs, 
having defined values of ric, Vc, l^c and lengths which lie 
in the small interval near Ic- In Eq. ( p7| ) we regard the 
statement c S P as a boolean function, which is equal to 
one if c g P and equal to zero if c ^ P. 

We can also count the number of pairs of the CPOs, 
c, c' e C with the fixed length difference Ic — Id , such 
that both of them are actually PPOs, c, c' E P. The 
probability of finding such a pair is defined as 



P(c,c') = (ce P n c' e P), 



P(c)P(c') 



i?2 i^c ; ) 



(28) 



This equation defines the averaging over ^"^2"' ■ '^^^ fluc- 
tuations of the Lyapunov exponent are ignorjCd here for 
two reasons: First, it was found numericallyt3 that the 
fluctuations of the stability amplitudes do not affect the 
correlation function; second, our result Eq. (|4^) is in- 
dependent of these fluctuations, because P(c) is smooth 
function of Xc. 

The probability that CPO c contains another CPO c' 
can be defined by averaging over the length of c 



P(c' G c) 



c e c 



^cpo(^c) 



This probability is related to the probability P{c) in 
Eq. (|^) since 



(30) 



where 1* is an arbitrary length such that lc/2 > 1* < Ic- 
Averaging over Ic gives 



+ E(-ir^' [nc &c)-\{p{pec) 



(31) 



where the composite orbits c' are such that 1^ > 1* and 
Vp e c' Ip < l*. We assume here that P(c' G c) = 
ripsc' ^iP S ^) unless 1^ is larger than I*. One can choose 

by minimizing the absolute value of the second term 
on the right hand side of Eq. (|3^) . 

The main purpose of this section is to construct an 
equation analogous to Eq. ( pl| ) for the joint probability 
P(c, c'). Equation (30) for pair of CPOs reads 



ce P n c' e P ^ f] p(^c n p^c' . (32) 

pGP, lp<l' 

This statement can be converted to a Boolean expression 
and averaged over {Ic + Ic')/'^- That is 



(33) 



^(c,c') = (np,z,</;[i-pec - pec' 
+ p e c n pec' 



Let CPO c" consist only of the PPOs which are shorter 
than^^. We can introduce the sum over such orbits into 



Pic,c') = (y2 n [i- pec -pec' 

c" p^c", lp<ll 

p£c" 



lr,<l* 



(34) 



^ E n [i - ^(P e c) - P{p e c') 

c" p^c" ,lp<l* 

X P(c" ec n c" ec') + Y[\i- Pip ec)- P{p e c') , 

(35) 



ip<it 



where the probability of having a "common divisor" is 
given by 

P(c" ecc^ c" e c) = U" e c r\ c" e c'\ . (36) 

We assumed that for p ^ p' 

(^pecn p'e c')^ « (p e c)^ (p' e c')^ ^ (37) 

and also we made the same assumptions as in the deriva- 
tion of Eq. (|3ll) . Particularly we neglected correction 
terms in Eq. (P5|), which are sums over c", Ic" > l*- Such 
correction terms were irnportant in Eq. (|3l|), but they 
can be neglected in Eq. ( |3q ) because the greatest contri- 
bution to the correlation of PPOs is given by the short 
CPOs c". 

It is convenient to rewrite Eq. (^5|) as a relation be- 
tween the correlation function of PPOs and CPOs having 
a "common divisor" 

P{c)P{c') -P{c,c') 

« E| [^(c" ^ '^)^(c" ^ c') - ^(c" e c n c" e c') 



n [i - G c) - p(p e c')] I 



(38) 



p<^c" , lp<l* 
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It is essential for the derivation of Eq. (38) that the com- 
posite orbit c" does not contain repetitions of the primi- 
tive orbits. This agrees with our definition of composite 
orbits Eq. (^) . In general one can allow repetition of the 
primitive orbits in the definition of the composite orbit. 
In this case Eq. (38) is also correct, but the sum must be 
taken over the composite orbits which does not contain 
repetitions of the primitive orbits. 

The derivation of Eq. (|3^) contains a number of ap- 
proximations. Therefore, one would like to verify that 
this equation preserves the normalization of the correla- 
tion functions. This verification appears in Appendix A. 

In certain cases, any given composite orbit contains 
a small number of PPOs and therefore we can make a 
further approximation in Eq. ( ^8|) 



n [i - Pip e - Pip e 

tc", lp<l' 

Jl [i - P(p e c) 1 - P{p e c' 



(39) 



lp<l'i 



The condition of a small number of "divisors" 



P{p £ c) (X e ^p'p < 1 

is not enough to justify making this approximation. It 
might happen that the error becomes larger when the 
sum over c" in Eq. (|3^) is computed. One can check that 
this error is of the same order as the contribution of the 
repetitions of PPOs which has been neglected. 

From Eqs. (|8|) and (|9|) by making use of Eq. ( ^ we 
obtain 

P{c)P{c') ~ P{c, c') « P(c)F(c') 

X ^ [p(c" e c)P{c" e c') - P(c" e c n c" e c')' . 

(40) 

The right hand side of Eq. (|8|) or Eq. (^) contains 
the non-trivial probability of finding a "common divisor" 
P{c" e c n c" e c'). it is difficult to find a "common 
divisor" if the lengths of c and c' are close to each other. 
We can compute this probability by noting that in or- 
der to find a "common divisor" of length Zc" one has to 
find two CPOs c"* and c™ of length l^ni = Ic — Ic" and 
l^ii: — Ic' — Ic"- The same is true for other components 
of Xc characterizing each CPO. The probability of find- 
ing a pair of composite orbits of specific lengths is given 
by the two-point correlation function of the composite 
length spectrum. Therefore, we substitute 



P(c" e c n c" e c') = P(c" e c)P(c" e c') 

--A ,,l n ^^{Xc — Xc" 1 Xc' — Xc 



(41) 



into Eq. (WG) and arrive at 



R2{Xc,Xc')- P{c)P{c') 
X ^e-^^"'="(72(A, 



Ac" , Xc 



Xc 



(42) 



In the more complicated cases, when the approximation 
Eq. fails, we should substitute Eq. (El into Eq. (Isf 



R2{Xc,Xc')^Pic)Pic') 

1 - P(p e c) - P(p e c') 



x< n T 



1 - P(p e c'] 



X 



E 



" C2 {Xc — Xc" , Xc' — Xc 



n 



pgc" 



1 - P(p e c) - P(p e c') 



(43) 



where the product in the braces is assumed to be con- 
vergent to some smooth function of both Ac and A^ but 
independent of the I* . 

The substitution of Eq. (p^) into Eq. (|2|) gives 



P2(A,A') 



(A)dppo(A') 
c'cpo(A)dcpo(A') Jq 



dx" 



J2 dcpo(A")C2(A-A",A'-A") 



(44) 



n" u" fi" 



where the summation in Eq. (^2|) is replaced by the inte- 
gration. 

Assuming that the product of the mean densities in 
Eq. (^) varies slowly on the scale of the A — A' de- 
pendence of the correlation functions we can apply the 
normalization conditions Eqs. (^ and ( p^ ) to get 



i^ppo(A) 



'^ppo(A) 



C?cpo(A) 



/•OO 

V / dx'dcpo(A')dcpo(A- A') 



(45) 



71' u' fl' 

which is satisfied by 



rfppo(A) cx 
dcpo(A) cx 7 



1 e 2m 



e 2x e 2a: 



y2Tix/a ^l-Kxjh 



(46) 



(47) 



The dependence on n is not known but assumed to be 
the same for both dppo(A) and dcpo(A). The parameters 
a and h here are lengths of the order of the billiard size. 
The mean density p£ PPOs Eq. (^6|) was introduced by 
Becj-y and Keating£3. The numeric tests of Dittrich et 
alEil support the distribution Eq. (p6|). Then the mean 
density Eq. ( |47| ) is_the solution of Eq. (^. 

Assuming that dppo(A)/dcpo(A) = (72;)"^ for large x, 
we obtain the main result of this section 



R(x,y) 



x^i^ Jo 



dx dcpo{x') C{x - x', y) , (48) 
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which was obtained by the summation on both sides of 
Eq. (^) over n, ^ at n' = n, v' = v, and /i' = /£• 

It is instructive to check how Eqs. (^) and ( pO| ) work 
for logarithms of integer and prime numbers. It turns 
out that Eq. (js^) is precisely equivalent to the Hardy 
- Littlewood conjecture. This can be shown by simple 
algebra and we put this calculation in Appendix B. At 
the same time Eq. (40 ) is equivalent to the smoothed form 
of the Hardy - Littlewood expression, see Appendix C. 



V. BEHAVIOR OF THE FORM-FACTOR NEAR 
HEISENBERG TIME. 



Ki;,{x,k)=-K\'{2nd{k)-x) 



diagV 



(52) 



This equation together with Eq. ( p9| ) expresses Koff{x, k) 
in the terms of ii'diagla;), -f^diag(^)' ^^'^ d{k). Equa- 
tions (|4^) and ( |5^ ) solve the problem of the periodic orbit 
computation of K{x^ k) near x ~ 27rd(/c). 

One can understand the important role of the short 
composite orbits by looking at another form of Eq. ( |49| ) 
which can be o btaine d fr om Eq. ( ^2|) by making use of 
Eqs. (|lOd|) and (p5^): 



Kosix, k) 



1 



2 ^ 



e-^^'^K^^s{x-l,,k) 



(53) 



For the case of billiards it is convenient to introduce 



the Heisenberg length 2'Kd{k). From Eqs. (10b) and ( p5q ) 
one can see that the diagonal parts of the form- factors are 
independent of the Heisenberg length, because they are 
independent of k. They are purely classical quantities. 

The information about the behavior of the form- factors 
near the Heisenberg length is hidden in their off-diagonal 
parts. The off-diagonal parts of the form-factors arc not 
independent; the Fourier transform of Eq. ( psf ) with re- 
spect to y gives us 



Kos{x, k) = 



1 



+ Kl.{x')Kl^{-x-x',k) 



Kijx')KUx^x',k) 

dx' , (49) 



wher e we substituted the mea n de nsity o f CP Os from 
Eq. (25c) and also used Eqs. (lOd) and (25e). In the 
derivation of Eq. ( ^ ) , we have assumed that the Heisen- 
berg length is so large that Eqs. (|lOd|) , (25c), ( ^5e| ), and 
( p8| ) are valid. We also put Eq. ([49|) in the form which 
is symmetric with respect to x. Equation ( ^ ) shows 
that Kos{x^ k) is an analytic function near the Heisen- 
berg length X ~ 27rd(fc) in agreement with results of Ref. 
p7| . A more detailed comparison of results is possible if 
we express K'^f^[x, k) in terms of -K^diag(a;). 

The functional equation for the spectral zeta function 
Eq. dl^ ) results in symmetry properties of the correlation 
function Eq. ([is|). The Fourier transform of this equation 
gives us 



(50) 



K'^{x, k) = K^{2TTd{k) - X, k) 



The derivative with respect to x of both sides of Eq. (|5' 
can be written as 



KiU^) + K;,{x) 



Kii^{27rd-x)-K'^^i27rd-x) , 

(51) 



where the argument k is omitted at d and at K^ff{x). 
Near the Heisenberg length K^^^^{x) is almost constant 
and its derivative can be neglected. The same thing is 
true for K^^{x) for small x. We, therefore, remain with 



valid for a; > 0. Substitution of the integrated Eq. (|5 
gives 



i^off(27rJ-a;) =^ 



47r7^ 



kLAx + Ic)~1 ■ (54) 



The function Kos{2TTd—x) has spikes for the small values 
of X equal to the length differences of the short CPOs. 

The inverse Fourier transform of Eq. (^9|) would give 
us the expression for i?off(e, k)- However, it does not con- 
verge, because K^ff{x, k) remains constant when x goes to 
infinity. The regularization of the inverse Fourier trans- 
forms obtained by taking derivatives of Eq. ( ^ ) with 
respect to x and 2'Kd gives 



i?off(e, k) 



cos{2TTd{k)e) 
2Tr^j^ 



|Cdiag(e)|'+rf(fc)5(e) . (55) 



This equation, which reproduces the results of Ref. |T^, 
was adopted for the ballistic systems in Ref. |3| and red- 
erived in Ref. ^ see Sec. ^ of present work for more 
detailed discussion. 

It is also instructive to check all of our equations for the 
case when the system exhibits universal behavior. The 
random matrix theory pfiedicted the correlation functiops. 
of the density of statesES and the spectral determinantES 
to be 



R{e, k) 
C(e,k) 



cos(27rd(fc)£) - 1 



27r2e2 



d(k)S{e) 



27 



Td(fc)e sin(7rd(fc)£) 



(56) 
(57) 



where the latter satisfies the symmetry relation Eq. (|is|). 
The commonly accepted form of Eq. ( |5^ ) can be obtained 
for R{e, k)/[d{k)]'^ and it is sm^{z)/z^, where z = TTed{k). 
The universal form-factors are 



K{x,k) 



2TTd{k) 



2TTd{k)\ 



87r2 



{x, k) = -fe{x)e{2TTd{k) - x) . 



(58) 
(59) 



One can note the symmetry of the universal form-factor 
Eq. (E^) with respect to the exchange x ^ 2iid{k). 



11 (November 30, 1997) 



In order to see some kind the of universal length spec- 
trum correlations, let us use the inverse density of states 
kx, which is a purely classical function defined in Eq. (|^). 
The Fourier transforms of Eqs. (^8|) and (|59| ) with respect 
to k and the change of sign give 



R{x,y) = dppo(a;)- 



dx 



I s m{k^,y) 
Try 



C{x,y) = dcpo{x) 



sm{kxy) 



Try 



(60) 
(61) 



Both correlation functions satisfy the normalization con- 
ditions Eqs. (0) and (p^). Equation ( [6l| ) is the "univer- 
sal" correlation function of the composite actions. Equa- 
tion (60) is the essence of the relation Eq. ( |4^ ) and it is 
the most general "universal" correlation function of the 
actions. 



VI. DISCUSSION AND SUMMARY 

The principal results of the present work are the sta- 
tistical relations Eq. (^) or E^. (|5^) and the approxi- 
mate functional equation Eq. { f)2\ } . The correlation func- 
tions and their form-factors depend parametrically on the 
widths of the averaging windows Afc and A/. At the same 
time these parameters do not appear in the statistical re- 
lations Eq. ( |49| ) and Eq. (|5^). Therefore, these relations 
are valid for some values of Ak and Al which should be 
specified. 

The probabilistic derivation of Eq. ( p9| ) means that 
Kos{x, k) and R{e, k) possess fiuctuations. Therefore, 
one has to choose Ak and Al in such a way that that 
these fiuctuations wilL be smaller than the correlation 
functions themselves H3 

Let us assume that the four-point correlation function 
of levels has the universal random mat rix t |heory form 
that is justified for some chaotic systemsoca. Then one 
can show, by making use of the geometric representation 
Fig. ^, that the mean square fluctuations of R{e, k) are 
of the order of J^/s A/fc-Vs. 

Let us assume that the two-point correlation function 
given by Eq. ( p5| ) behaves like for ed{k) 1 . Then 
it becomes of the order of AP near e ^ 1/AZ, which 
is the maximal value of e where the definition Eq. (^ is 
meaningful. Therefore the condition that the fiuctuations 
of ii(e, k) are smaller than R{£, k) itself is 



AP > 



Afci/2 ■ 



(62) 



This inequality, together with Eqs. (|3|) and (|lj), re- 
stricts the choice of Ak and AZ strongly. 

The inequalities Eq. (|l^ taken near the Heisenberg 
length become the approximate equation relating Al to 
Afc. We can substitute A/(Afc) into Eq. (p2) and obtain 



A, 5 



(63) 



Let Af{k) be a function counting energy levels with fc„ < 
k. One can apply the approximation d'{k) « (P{k)/Af{k) 
which is valid for generic system. The counting function 
can be substituted into Eq. ( |63|) , and it gives the widths 
of the averaging window measured in the number of en- 
ergy levels 



AM » . 



(64) 



This inequality has to be satisfied together with AJV <C 
Af. Therefore the numeric check of relations like Eqs. (m) 
and (^5|) is difficult. One has to take at least A/" ~ 



and average over AJ\f 10^ levels in order to see non- 
universal features predicted by Eqs. ( |5^ and (^5|). The 
two-point correlation function has to be evaluated for e 
in the range from zero to d/Al ^ Af/AJV ^ 10 level 
spacings. This estimate also shows that 



AkAl 



AAf^ 



> 1 



(65) 



and the definitions of the correlation functions Eqs. @) 
and ( p^ ) are justified. 

The basic object characterizing a chaotic dynamical 
system is the set of mixing rates. They show how fast 
the density-density correlations decay. These rates are 
zeros of the so-called classical zeta function, see e.g. Ref 
|36| , which is approximately Cdiag(is), see Ref. The 
functional equation in the form Eq. (b3) implies 



eCoff (£, fc) = e-^--^^)^ [£Cdiag(£)]' 



(66) 



where eCdiag/off is defined as the inverse Fourier trans- 
form of ^'di'ag/off- Therefore Coff(£, k) is the classical 
zeta-function modulated by quantum oscillations. This 
function contains the information not only on the mixing 
properties of the billiard, but also on quantum recur- 
rence. 

The association of Cdiag(is) with the classical or Ru- 
elle zeta function is valid in the same approximation as 
Eq. ( ^ was derived. Therefore, Eq. (|5^) reproduces ex- 
actly results of Ref. for the case of the system with 
broken time reversal symmetry. The result of Ref. ^ goes 
beyond the range of lengths given by Eq. (|lT|). One has 
to make the replacement in Eq. (pq) 



l^?diag(£)|' 



[27rd(fc)7]^ 



C(fc + £/2)C(fc-£/2) 
C(fc + e/2)C(fc-e/2) 



diag 

(67) 

and use the product Eq. (|l7| ) truncated near the Heisen- 
berg length Ip < 2Trd(k) in order to obtain the more ac- 
curate results of Ref. 0. We performed the more accurate 
computation for the prime numbers, see Appendix B, and 
obtained the precise form of the Hardy - Littlewood con- 
jecture, which is equivalent to the replacement Eq. ( |6^ ) 
in Eq. (||). 

All calculations in this work were done for chaotic bil- 
liards, where the action of orbit p has the very simple 
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form Sp = hklp. The present theory can be general- 
ized for other types of dynamical chaotic systems. The 
only problem is that the correlation function of actions 
(^6{x — Sp)S{x' — Sp')) has to be defined by averaging 
over the constant mean period dSp/dE + dSp'/dE. The 
^stification of such procedure was discussed in Refs. ||, 

The result of the present work cannot be applied to 
the systems which possess time reversal symmetry. In 
this case we should compute C2{X',^') or R2{X;X') 
for n ^ fi' , i.e. the action correlations of orbits with 
different winding numbers. However, the results pre- 
sented here can be generalized for a system having dis- 
crete symmetriestd, see also a discussion of-this problem 
by Leyvraz and SeligmanEZI and Agam et aEB. In the case 
of the discrete symmetry of the system, one should com- 
pute the multiplicities of the periodic orbits and correla- 
tions of these multiplicities in order to obtain the spectral 
statistics. Probabilistic number theory methods can be 
used again for this purpose as in the case of the modular 
group. LJ 

In summary we have shown that the functional equa- 
tion for the spectral determinant implies the action corre- 
lations. Quantum-classical time scale separation for sys- 
tems with broken time reversal symmetry allows us to 
compute all spectral and action correlation functions be- 
yond their universal random-matrix theory shapes. One 
of the central technical points in the present work is the 
derivation of the Hardy - Littlewood conjecture of prime- 
number correlations in such a way that it can be used for 
actions of a dynamical system. 
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APPENDIX A: NORMALIZATION OF THE 
CORRELATION FUNCTIONS. 

The normalization conditions Eqs. (Q) and ( p^ ) are 
generic for the correlation function and we expect that 

E,.[P(c)P(c') - P(c,c')] - Pic) and Ec' [^(c" ^ 

c)P(c" e c') - P{c" e c n c" e c')j = P{c" e c). The 

sum over c' in Eq. pq) leads to 



We can represent P{c" G c) as a product, and then 



(Al) 



■ ■ c" pec"jp<i'^ ^-^ ' 

Pip e c) 



i„<i' 



n [i - 2Pip e c) 
n [i - Pip e c) 



1 



1 - 2P(p e c) 



in agreement with Eq. (|3l| 



APPENDIX B: CORRELATIONS OF PRIMES 
OBTAINED FROM CORRELATIONS OF 
INTEGERS. 

We can demonstrate the relation between the correla- 
tion functions of composite and primitive actions using 
the example of integer and prime numbers. Let us con- 
sider the prime number p as a PPO of length log(p) and 
the integer number n as a CPO of length log(n). This 
"CPO" may contain repetitions of "PPO" , but this is not 
important for large numbers. 

Saying that the "CPO" m "is a part of" "CPO" n 
means that to is a divisor of n, and we will write m\n. 
This notation allows us to rewrite Eq. (Bq) as 



P(n)P(n') - P(n,n') 

= ^ j |^P(to I n)P(TO I n') - P(to \n r^ ■m\n') 



n [l-Pip\n)-Pip\n') 



(Bl) 



where m runs over integers which do not contain powers 
of primes. Here P(ri) is the probability to find a prime 
number, and Pin, n') is the probability to find pair of 
primes, see Refs. mi M for precise definitions. 

The derivation of the Hardy - Littlewood expression 
is based on the probabilityCj that m is a divisor of n 



Pirn I n) 



(B2) 



Our derivation is based on the probability to find a com- 
mon divisor 



P(m I n n TO I n') = — (5- 

m. ^ — ^ 



n~n' —ml j 



(B3) 



where the Kroneker 6 symbol in the right hand side is 
equal to one for n = n' — ml and equal to zero otherwise. 
It carries information about correlations of integers, be- 
cause if we found with probability 1 /to that n is multiple 
of TO, then all integers n' = n + ml are also multiples of 
TO with the probability 1. 

Substituting the probabihty Eq. ( p^ ) into Eq. (Bl) we 
obtain 
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P(n)P(n') - P{n,n') 

-i:{[;^-;!;E*..,..-.„in[i-^l}- (B4) 



and the right hand side of this equation is zero if m is odd. 
Indeed, if m is odd then p = 2 is not a divisor of m and 
1 — 2/p gives zero. The density of prime numbers, P(ji) = 
Y\p{l — must foUow the prime number theorem and 
this is the criterion for choice of the upper hmit in the 
products over primes. 

The next step, which is not necessary, is to separate 
the smooth part from the part containing correlations 



P(n,n') = 



E 



1-^ 

m ^ L p. 

m I (n— n') pf rn 



(B5) 



and this expression is non-zero only if n— n' is even. Only 
terms with even m contribute and the term with m = 2 
is always present. We can, therefore, extract the factor 
2 from all terms, and we obtain for even n — n' 



p>3 

+ E 

m I (n— n') 



2 

1 - - 
P 



1 
2m 



n [> 



(B6) 



where m is odd and does not contain powers of primes. 
We can rewrite the sum as a product 



p>3 



P 



n 



p I (n— n'), p>3 



p-i 

p-2 



(B7) 



This product contains the correct enhancement^ of the 
probability to find a pair of primes n,n' by the factor 
{p — l)/{p— 2) per each prime, which is divisor of n — n' . 
A similar computation shows that 

m vl'm P 



Therefore we have shown that Eq. (B5) follows Eq. (B4). 

The product Eq. (B7) diverges to zero for large n. For 
this reason the probability P{n, n') is divided by P{n) — 
— 1/p), seeRef.^, Chap. Postscript on prime-pairs. 



P{n,n') 
P{n)P{n') 



l-2/p 



p>3 



(1 - 1/py 



n [i 

p I (n— n') 



i/p 

\~2Ip 



(B9) 



This correlation function contains only the convergent 
products and it allows one to compute statistics o£-.the 
zeros of Riemann zeta function with great accuracyQ 



APPENDIX C: SMOOTHED CORRELATION 
FUNCTION OF PRIMES. 

The smoothed form of the correlation function of prime 
numbers is valid for \n — n'| ^ 1, and it can be obtained 
directly from Eq. 



P{n,n') 
P{n)P{n') 



— 2 ^ P{m \ n C^ m \ n') 



(CI) 



In this equation we used the probabilities of finding the 
prime number P{n) and the probability of finding the 
pair of prime numbers P{n, n'), see Refs. and we 

substituted instead of P(m | n)P{m\n'). 

The smoothing of the probability to have a common 
divisor can be done in a number of different ways. We 
suggest that the simplest one is 

1 /""/^ 1 m 

— / P(m I n n m I n')dn' = -^0{\n - n'\ ) . 

m J-m/2 2 

(C2) 

|) can be replaced by 
l-e(\n-n'\ - to/2) 



The summation over to in Eq. 
the integration and we obtain 



P{n,n') 
P{n)P{n') 



dm ■ 



1 



2 71 - n' 



(C3) 



The choice of the lower limit in the integral is not impor- 
tant, because the integration is performed from 2|n — 
This result is the leading order expansion of the correla- 
tion function in |n — 7i'|^^ and it coincides with Keating's 
result EJ 

The probability of finding a prime number can be com- 
puted from Eq. ( |C3| ) by application of the normalization 
condition Eq. 



1 



P{n) 



E 



P{n,n') 
P{n)P{n') 

1 _ 

2\n-n'\ ^ 



log(n) 



(C4) 



\n—n' I <n 

according to the prime number theorem@@. 
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